Bounce of gravitons emitted from the brane to the bulk 
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Abstract 



' Solutions of geodesic equations describing propagation of gravitons in the bulk are 

studied in a cosmological model with one extra dimension. Brane with matter is embedded 

in the bulk. It is shown that in the period of early cosmology gravitons emitted from the 

brane to the bulk under certain conditions can return back to the brane. The model 

is discussed in two alternative approaches: (i) brane with static metric moving in the 

AdS space, and (ii) brane located at a fixed position in extra dimension with non-static 

| metric. Transformation of coordinates from the one picture to another is performed. In 

both approaches conditions for gravitons emitted to the bulk to come back to the brane 

C*") . are found. 

> 
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^ ! 1 Introduction 



Models of the space-time with extra dimensions beyond the four macroscopic can be considered 
as candidates for alternative cosmologies. The minimal requirement on these models is an 
ability to reproduce the main observational cosmological data such as the age of the Universe, 
abundance of light elements produced in primordial nucleosynthesis, etc. 

In this paper we study a 5D model with one infinite extra dimension with the 3-brane 
with matter embedded in the bulk [1, 2, 3]. In the radiation-dominated period of the Universe 
interacting particles of hot matter on the brane produce graviton emission which can escape to 
the bulk. This effect alters time evolution of matter on the brane as compared to the Standard 
cosmological model, and, in particular, alters primordial nucleosynthesis [4, 5, 6, 7]. We study 
solutions of geodesic equations for gravitons propagating in the bulk. In the period of early 
cosmology, in which in the generalized Friedmann equation the term quadratic in brane energy 
density is dominant, there appears a new effect: gravitons emitted in the bulk can bounce back 
to the brane. 
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In a class of metrics defined below, the 5D model can be treated in two alternative ap- 
proaches. In the first approach the brane is located at a fixed position in the extra dimension, 
and the metric is time-dependent. In the second approach the brane is moving in the static 5D 
space-time, constructed by attaching two AdS spaces at both sides of the brane. We consider 
a method to connect both approaches (cf. [8, 9, 10]). 

Solving the geodesic equations in the period of early cosmology, in both approaches, we 
find solutions for trajectories of gravitons propagating in the bulk and conditions for emitted 
gravitons to bounce back to the brane. Related issues were discussed in [11, 12, 13]. 



2 Static and time-dependent netrics 

In the first approach we consider a class of non-static metrics [1, 2] 

ds 2 = -n 2 (y, t)dt 2 + a 2 (y, t)dx a dx a + dy 2 . (1) 

The components of the metric n 2 (y,t) and a 2 (y,t) satisfy the system of 5D Einstein equa- 
tions with the energy-momentum tensor of matter confined to the brane and the cosmological 
constant A. The brane is located at a fixed position in the extra coordinate, which we set 
y = 0. The function n(y,t) is normalized by condition n(0,t) = 1. For simplicity we consider 
a spatially flat brane. 

The energy-momentum tensor on the brane is 

T v = diag{-p-a,p-a,p-a,p-a}, (2) 

where p and p are the energy and momentum densities of matter on the brane, a is the tension 
of the brane. 

It is convenient to define the normalized energy density, momentum and tension of the brane 

2 " 2 " 2 " 

K p K If K a 

where p and a have dimension [mass] 4 and p and a have dimension [mass] 1 . Parameter p 
having dimension [mass] is defined via the 5D cosmological constant p 2 = k 2 \A\/6, where 
k 2 = 8it/M 3 is the 5D gravitational coupling constant. 

In the class of metrics (1) there is a unique extension of the metric from the brane to the 
bulk [2] 



2, ^ « 2 (0,t) 



a 2 (y,t) 



e 




(3) 



and n(y,t) = a(y,t)/d(0,t). The function a(0,t) satisfies the generalized Friedmann equation 
(7) [1, 2] 

H 2 {t) = -p 2 + (p + a) 2 + pp w {t), (4) 
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where 

H(t) = ^4, Pw (t) Pw0 



a(0,ty ,wy ' a 4 (0,t)' 

are the Hubble function and the Weyl radiation term. In derivation of (4) p w0 appears as an 
integration constant [2]. 

In the second approach the moving brane separates two static 5D AdS spaces attached at 
both sides of the brane with the metrics [14, 15, 16, 17] 

ds 2 = -f(R)dT 2 + — — + p 2 R 2 dx a dx a = gijdx'dx*. (5) 

Here 

f(R) = p 2 R 2 - ^. 

R and T have dimension [mass]^ 1 . Reduction of the metric on the brane is ds 2 = —dt 2 + 
p 2 R 2 (t)dx a dx a , where t is the proper time on the brane. 

In a parametric form trajectory of the brane is defined through the proper time t on the 
brane as R = R b (t), T = T b (t), x a = 0, where -f(R b )f 2 + f-^R^R 2 = -1. It follows that 

where £ = ±1. 

From the spatial components (a,b — 1,2, 3) of the junction conditions on the brane follows 
the generalized Friedmann equation [1, 2, 14, 15, 16, 18, 3] 

f) =V + (p + -) 2 + ^. (7) 

Kb I K b 

Equation (7) is of the same form as (4). The terms a 2 (0,t) and p w (t) = p W Q/a A {t) can be 
identified with the terms p 2 R 2 (t) and P/pR^it). Below we set a — p and consider the case 
with P = [7]. 



3 Connection between two pictures 

To establish connection between two pictures, we consider geodesic equations. In the picture 
with the static metric, geodesic equations with parameter y along the geodesies are 

d 2 T rp dTdR 
d 2 x a ^ dx b dR 

^ + r U^)\r ?T ( d S)\r^ d 4^, ( io> 



dy 2 \dy J \dy) a0 dy dy 
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where the Christoffel symbols are 

V T _ J V R _ J V R _ 1 f f t V R _ f,,2p -pa _ °b 

~ 2/' R ~ ~2f' 1 T ~ 2 a6 ~~ 6 _ i?' 

Here (T, i?) = (T ± , i? 1 * 1 ) are coordinates in the AdS spaces at the opposite sides of the brane 
/ = fi 2 R 2 . Integrating the geodesic equations, one obtains 

<IT ' E " (b '" C " f(W±V = f(R)(C R ±) 2 + E± 2 -C« 2 , (11) 



dy /(-R)' dy \i 2 R 2 ' \ dy 

where (E ± , C a , C R± ) are integration parameters. (dT/dy, dR/dy, dx a jdy) are the components 
of the tangent vectors to geodesies which we normalize to unity. Imposing the normalization 
condition 

dx % dxi 

1 — J-9ij = l, i = T,R,a 
dy dy 

we obtain that (C R± ) 2 = 1. 

Connection between the two pictures is established in the following way. We consider 
solutions of the geodesic equations even in y: T + (y) = T~(—y), R + (y) = R~(—y) . 

The hypersurface (T, R, x a = 0) is foliated by geodesies that intersect trajectory of the brane 
(Tfe(t), Rb{t), 0) and at the intersection point are orthogonal to it (cf [9]). Geodesies are labeled 
by parameter t. t is constant along a geodesic. The geodesies orthogonal to the trajectory of 
the brane satisfy boundary conditions at y = : R ± (0,t) = Rb(t), T ± (0,t) = Tb(t). Setting 
C a = 0, have 

0r*(y,t) Ee(y) dR±(y,t) - 2] i/2 

where a — ±1, and E + = E~ = E. 

The normalized velocity vector of the brane and the normal vector to the brane trajectory 

are 

4 = (f b , R b ) = ^f(R b ) + R 2 /f(R b ), R b )^J = V e(y) (^j^- y ^ f(R b ) + R 2 

'(13) 

where r\ = ±1. From (12) we obtain the tangent vector to a geodesic 2 

ut± = ( W «<v)y/fW + &) ■ ( 14 ) 

By construction, the tangent vector to a geodesic at the intersection point with trajectory of 
the brane is (anti) parallel to the normal to the trajectory of the brane, 

u l | w=0 | K (15) 

From this condition it follows that 

E = V R b , a = £r] (16) 



2 On the brane we find n R = (dR/dy)n v and n T — (dT/dy)n v , where n v is the normal to the brane in the 
metric (1). The same relations hold for v\. 
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4 Integration of geodesic equations 

Let us integrate the geodesic equations. From the second equation (12) we obtain 



R 2 

R ± (y,t) = R b (t) cosh/iy + a\j — | + R 2 sinh/x|y|. 



Using the Friedmann equation (7) with P = 

H 2 = Rl/Rl =p 2 + <l W 
and introducing j3 = sign(R b (t)), we rewrite (17) as 

R ± (y, t) = Rb(t) ^cosh fxy + af3 ^1 + — ) sinh ji\y 

For correspondence with (3) (the case with P = 0) 



a 2 (i/,t) 



a 2 (0,t) 



; 2^| M 2 + e -2,\y\ (P +2 ] 2 _ 2 P(P +2 

\fij \n J n\n 



we set 



ad = -I. 



Omitting in (17) (±), we obtain 

R(y,t) = 

Introducing y , such that 

we express R(y, t) and R'(y, t) as 
R b {t) 



R b (t) 





_ e n\v\P 




P. 



P 



p + 2/i 



1/2 



R(y,t) 



Integrating Eq. (12) for T(y,t) with R(y,t) (20), we obtain 

r± 1 coshWj^ + 

fiE smh(n\y\ + fj,y ) 



(17) 



(18) 



(19) 



(20) 



(21) 



sinh(/x|j/| + /xy ), t) = -e(y)R b (t) cosh(/x|y| + /xy ) (22) 



(23) 



For C + (t) = C (t) — C(t) the limits y = from both sides of the trajectory of the brane 
coincide. 
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To determine C(t), first, we consider transformation of the metric (5) with P = from 
coordinates R, T to y, t. We have 

ds 2 = dy 2 (-p 2 R 2 T' 2 + -^-] +2dtdy { —p 2 R 2 TT' + ) +dt 2 { -p 2 R 2 f 2 + ^— \ +p 2 R 2 dx a2 . 



p 2 r 2 J * y ^ p 2 R 2 J \ v 2 R 2 

(24) 

Substituting T' and R' (12), we find that the coefficient at dy 2 is e 2 (y). The coefficient at dydt 
is zero, if 

• = RR' = RR' 

p 4 R 4 T' ii 2 R 2 {t{y)EY 1 ; 

where in the second equality we have used (12). Substituting in the coefficient at dt 2 expression 
for f (25) and using (12) and (22), we obtain R 2 / R b = n 2 (y,t) (cf. (6)). 
Writing (23) as 

and taking the time derivative, we obtain 

f- kR ' 1 d ( R <y)\ | g (27) 

p 2 R 2 (e(y)E) p 2 R dt\ E J K ' 

The first term in the rhs of (27) is the same as in (25). Remarkably, substituting explicit 
expressions (22) for R and R', we find that the second term in the rhs of (27) is independent 

1 d (R'e{y)\ _r)y (t) (2g) 



p 2 Rdt\ E J R b (t) 

Choosing 



C = nf, (29) 
Kb 

we obtain T as in (25). In the limit y = we have 3 

T(y,t)\ y ^ = -r,£±£ = f b (t). (30) 
fi K b 

In the radiation-dominated period, conservation equation for the energy density on the brane 
is 

p = -4pR b /R b . 

Using this equation, we have 



4p • p + p 



3 Conversely, requiring that the limit y — of (25) is Tb(t), we arrive at (29). 
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From (20), using the conservation equation for p, we obtain 



R(y,t) 



Rb 



e"" 1 - e 

p 



(31) 



Jacobian of transformation from T, R to t, y is 



J = TR' - RT' = er]S-. 

Rb 



Although both R and Rb are non-zero, the expressions for T and T' contain R(y, t) in denomi- 
nator. R(y,t) is zero for e 2 ^' — (p + 2p)/p, and transformation from (T, i?) to (t, y) is valid 
for e 2 ^l < 1 + 2p/p. In the region e 2 ^ 1 > 1 + 2p/p R{y,t) can be defined as 



R(y,t) = 



R b (t) 



p 



*\v\ (E + 2 
p 



which is positive and increasing function y. 



5 Null geodesies 

Next, we consider null geodesies along which gravitons propagate in the bulk. In the picture 
with the static space-time null geodesies again are found from Eqs. (8)-(10). The tangent 
vectors to a null geodesic satisfy the relation 

dx % dxi , . 

* dX ix = °' (32) 

where, to distinguish the case of null geodesies, we relabeled the affine parameter from y to A. 
Substituting solution (11) in (32), we obtain that C R = 0. Taking for definiteness direction 
A > 0, we obtain tangent vectors to a null geodesic 

dT C T dx a C a dR T r- 

dx^JiRy H = ^ lx =ulc lv ^^' (33) 

where v — ±1. Here 

7 = C a2 /C T \ 

The geodesic equations can be integrated in a closed form. 

C T A 



R(X, t) - R(0, t) = y^l-^A, T(A, t) - T(0, t) 



(34) 

Here t is proper time on the brane at which begins the geodesic. 
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Let (-R(A), T(A), x a (X)) is a point on a null geodesic. R, T can be expressed as functions of 
y, t on the hyperplane x a = 0. For a point on a geodesic we have R(X) = R(y(X),t(X)), T(A) = 
T(y(X),t(X)). 




Figure 1: Schematic plot of geodesies on the brane and in the bulk. Curve 1 is the intersection 

of the brane world scheet 5 with the hypersurface R, T (brane trajectory). Curves 2 represent 

the geodesies orthogonal to brane trajectory in the hypersurface R, T. At the brane trajectory 

y = 0, and t is running along trajectory. At the geodesies orthogonal to the brane trajectory 

t = const and y varies along geodesies. Curve 3 is a null geodesic in the bulk bouncing back to 

the brane. Curve 4 is projection of the null geodesic on the hypersurface R, T. 

Let us verify that the vector (dy(X)/dX, dt(X)/dX, dx a /dX) is null in the time-dependent 

metric (1) 



where a 2 (y,t) = fi 2 R 2 (y,t) and n 2 (y,t) = 
From the equations 

dT(X) 

dX 
dR{X) 



tx) -"^(JD +a2 ^ 

R 2 (y,t)/R 2 (t). 



dX 



dT 
~dt 
dR 
~dt 



dt dT 
dX dy 



dt 

AX + 



dR 
dy 



dx a 
~dX 



dy_ 
dX' 

dy 
. dX 



0. 



(35) 



(36) 



we obtain 



dy _ RdT/dX-f dR/dX 
dX~ RT'-TR' 
dt _ T'dR/dX- R'dT/dX 
dX~ RT'-TR' 



(37) 



Substituting expressions (37) in (35), we rewrite this relation as 



f)rii+(£)>-i* 2 > 

2 2 

dx a \ ? 



^2^^(Rf-S-rRi) + p 2 R 2 (RT'-fR') 2 ( 
dXdX \ R 2 V 



Using the expressions (25) for f, (20) - (??) for R, R', (12) for V and (33) for (dR/dX, df/dX, dx a /d\), 
we verify that the relation (38) is valid. 

In the setting with the non-static metric (1), the system of the geodesic equations is com- 
plicated, and it is difficult to integrate it in a closed form. However, the equation for x a 

d 2 x a dx b dy dx b dt 

IX 2 + ay ~dX dX + bt llXdX 

can be easily integrated and yields 

dx a C a 



dX a 2 (y,t) 

which is the same as in (33). 



6 Conditions for bounce of gravitons 

Let us solve the combined system of equations for brane trajectory and a null geodesic and 
show that under certain conditions null geodesic can bounce back to the brane. We consider 
the radiation-dominated period of the early cosmology, when p > p. We assume that solving 
Friedmann equation we can neglect the energy flow from the brane to the bulk 4 . We obtain 

Apt + 8 (/it) 2 , (39) 



tit) 



and H/p = 1 + fit. 

As in [5], the equation for the brane trajectory can be written as 



where (5 = sign(Rb(t)). Substituting H 2 = p 2 + 2pp and expanding the rhs of (40) in powers 
of p/p < 1, we obtain 



4 As it is seen from (40), it is sufficient to assume that the terms containing bulk energy- momentum tensor are 
small. Expanding expression (40) in powers (i/p, we obtain that in the leading order these terms are canceled. 
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Using that p(t) = p(ti)(R h (ti) / R b (t)Y , we have 



2 P y - 2 P 2(t,)^(( 1 )- 

Integrating Eq. (41) with the initial conditions Rb = Rb(ti), T5 = Tf,(ti), we have 

evpiw - 7M*,)) = + ^ + j^Smr {m ~ Rl(u)) ■ (42) 

Here p? j p 2 (t\) is a small parameter. 

Let us consider a graviton trajectory in the picture with the static metric. From the equa- 
tions (33) we obtain 

§ = ™(l- 7 )Vy^ (43) 
al 

where k = sign{C T ). Integrating this equation with the initial conditions R = Rb(ti), T = 
T b (ti), we have 

TTTTT " i = ™(1 - 7) V V(T - Tift)). (44) 
^WlJ it 

To find, if at a time t the graviton trajectory returns to the brane world sheet, i.e. R = Rb{t) 
and T = T b (t), we equate solutions (42) and (44). In the case £/3 = z//t we have 

Hi - 7) -v» _ n (_J i-^j = _Jf! L_ f (ML) 7 _ A . (45 ) 

Eq. (45) means that graviton emitted from the brane at time t\ has bounced back on the brane 
at time t. 

Using relations a(3 = —1 and a = £77, condition £/3 = z//t is expressed as 

urjK, — — 1. (46) 
Setting z = Rbit^/Rbit) < 1 and 7/2 = [(1 - 7 )~ 1/2 - 1], we transform (45) as 

(1 - J) _ 7 7P 2 (t 1 ) 



z 7 {l-z) ' p 



2 



(47) 



The function of z in the l.h.s. of (47) is monotone decreasing with the minimum at z — 1 equal 
to 7. Eq. (47) has unique solution, provided 7 p 2 (ti)//i 2 > 1. Expressing p(ti) through p(t), we 
obtain 

^(l-^ 7 ) 7 7P 2 (t) 

"T=T~ " (48) 

The function of z in the l.h.s. of (48) is monotone increasing with the maximum at z — 1 equal 
to 7. Eq. (48) has unique solution provided 

ip 2 (t)/p 2 < 1. (49) 
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Existence of solution means that graviton emitted from the brane at proper time t\ comes back 
to the brane at time t. 

Momentum of graviton is proportional to tangent vector to a geodesic (33). In the basis 
(er» e «> e «)' where 

4 = (1/^,0,0,0,0), 4 = (0,^,0,0,0), et = (0,0,l a /fiR) 

the lengths of momentum components are 

|/| 2 = \g TT p T2 \ = C T2 /n 2 R 2 , \p R \ 2 = C R2 /p 2 R 2 , \p a \ 2 = C a2 /p 2 R 2 . 

In this basis condition (49) means that \p T \ ~ \p R \ 3> \p a \ 5 . 



6.1 Bounce of gravitons in time-dependent metric 



Let us show how the bounce of gravitons takes place in the picture with the metric (1). The 
difficulty is that in this metric we have no analytic solution of geodesic equations. We obtain 
equation for graviton trajectory using solution of geodesic equations (33) and transforming it 
to coordinates (t, y). 

From the equations (37) we obtain the equation of trajectory 



dy 
dt 



RdT/dX-f dR/dX 



V dR/dX- R'dT/dX 
Using (33), denominator in (50) is obtained as 



(50) 



T'dR/dX- R'dT/dX = 



C T R b e{y) 
p 2 R 2 



ur]KH(l - 7 ) 1/2 + 



(p + 2fi)e'^ + pe 



Ml 2/ 1 



It is seen that the expression in brackets is positive. 
The numerator can be written as 



RdT/dX-fdR/dX = 



Eq. (50) takes the form 



dz 
jidt 



1 



,3p 



C T R 
fi 2 R 2 



3p 



1 + vt]k{1 — 7) 



1/2 



e -»\y\(p + 2ii) + e^p 



2H 



(51) 



2H + VT^/T^j^jp + 2/i) + zp] 
2Hy/T^ + vtik\z- x (p + 2p) + zp] ' 



(52) 



where z = e^l In correspondence with (46) we take vr\ k — — 1. 
The numerator has zeroes at the points 



*l,2(f) 



1 + 



2p 
P 



1/2 



i_±vr\ 1/2 



(53) 



3 In different bases condition (49) written in terms of components of momentum may take different form. 
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and is positive for z 2 < z < z 1 . At the initial time t 1 y(t\) = and z(ti) = 1. Condition that 
z 2 {t) «g 1 is 

p 



/i + p(t) > 

The above formulas are valid in the region where R(y, t) > 

i<«,.(i + af. 



Eq. (52) can be written as 



3p 
p 



- 2 



(z- Zl )(z-z, 



Introducing new variable u = z — z 2 , we have 
"3p 



/i 



((z 2 + u) 2 - 1) + 2) 



(z - Zl )(z - z 2 ) + 2H 1 z/p y /T^ 
— (zi — z 2 — u)u 



7- 



(54) 



- z 2 - u) + 2^z(z 2 + -u)/a/1 - 7 



; (55) 

Expanding all the terms in (55) to the linear order in u and in small terms ypy and p/p, we 
obtain 



du 
pdt 



3p ( 2p 
p 



-2^ + 2u +2 



"fr- 2 Vfl 4 fl-^-^V. (56) 



At the initial time t\ at which trajectory begins, ypy > p/p and ^(^i) < 1- At larger times the 
point z 2 (t) eventually moves to the right and crosses unity at time t 2 . At times t > t 2 we have 
y/^y < fij p and z 2 (t) > 1. The initial time ti at which trajectory begins can be chosen so that 
22(^1) is arbitrary close to 1: 1 — z 2 (ti) — u(ti) <C ^/t". In this case £ 2 — ^i is also small and 
u(t 2 ) < v/7. 

Eq. (56) is simplified as 



du 
fidt 



^ f^-2^ + 2^+2 
fi \ p 



u 



4 = 



3p + «(8/i/3p - ^7) - 4/1^/3^ 



As a function of u the rhs of (57) is negative in the interval (u±, u 2 ) and y/j > u, 



(57) 



«1,2 = g 



Here «i > and u 2 < 0. 

Eq. (57) can be written as 



du 
fidt 



3p 2u-^y 
— u — 4- 



fX 



yfi-U 
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For u <C y/j, substituting p/p ~ Apt, this equation can be simplified as 



du 3u 
pdt Apt 



-4. 



(58) 



In the case ^fy 3> u solution is 




3/4 



16 ft 



7/4' 



u(t) 



u(t) 



(59) 



where t is a reference time after which approximation is valid. It is seen that u(t) — is reached 
in finite time t. For larger times t > t from (54) it follows that dz/ pdt < 0, i.e. the distance 
z(t) from the brane decreases. 

To conclude, in two approaches we have shown that under certain conditions emitted gravi- 
tons moving in the bulk along null geodesies can bounce back to the brane. 

I thank my colleagues at the Skobeltsyn Institute for discussion. My special thanks are to 
I. Tuitin and B. Voronov and M. Smolyakov for useful remarks. 
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